Relativistic Hartree-Fock-Bogoliubov model for deformed nuclei 



O 

(N 



(N 



o 

=5 



(N 
> 

o 

(N 

d 



J. -P. Ebran, 1 E. Khan, 1 D. Peha Arteaga, 1 and D. Vretenar 2 

1 Institut de Physique Nucleaire, Universite Paris-Sud, IN2P3-CNRS, F-91406 Orsay Cedex, France 
2 Physics Department, Faculty of Science, University of Zagreb, 10000 Zagreb, Croatia 

(Dated: July 28, 2011) 

The Relativistic Hartree-Fock-Bogoliubov model for axially deformed nuclei (RHFBz) is intro- 
duced. The model is based on an effective Lagrangian with density-dependent meson-nucleon cou- 
plings in the particle-hole channel, and the central part of the Gogny force is used in the pairing 
channel. The RHFBz quasiparticle equations are solved by expansion in the basis of a deformed 
harmonic oscillator. Illustrative RHFBz calculations are performed for Carbon, Neon and Magne- 
sium isotopes. The effect of explicitly including the pion field is investigated for binding energies, 
deformation parameters, and charge radii and has an impact on the nuclei's shape. 

PACS numbers: 24.10.Cn, 24.10.Jv, 21.60.Jz, 21.30.Fe, ll.10.Ef, 21.10.Dr 



I. INTRODUCTION 



Among the microscopic approaches to the nuclear 
many-body problem, nuclear energy density functionals 
(EDF) represent a tool of choice for the description of 
both static and dynamic properties of nuclei over the 
whole nuclide chart. EDFs subsum nucleonic short-range 
in-medium correlations, whereas static long-range corre- 
lations (deformation, pairing, ...) are incorporated by 
allowing a single-determinant state to break the sym- 
metries of the nuclear Hamiltonian [1]. A variety of 
structure phenomena in stable and exotic nuclei have 
successfully been described by EDFs based on the non- 
relativistic Gogny and Skyrme Q effective interactions, 
as well as on relativistic phenomenological Lagrangian 
densities Q. There are significant advantages in using 
covariant functionals [||. The most obvious is the natu- 
ral inclusion of the nucleon spin degree of freedom, and 
the resulting nuclear spin-orbit potential which emerges 
automatically with the empirical strength in a covariant 
formulation. The consistent treatment of large, isoscalar, 
Lorentz scalar and vector self-energies provides a unique 
parametrization of time-odd components of the nuclear 
mean-field [f|, i.e. nucleon currents, which is absent in 
the non-relativistic representation of the energy density 
functional. The empirical pseudospin symmetry in nu- 
clear spectroscopy finds a natural explanation in terms 
of relativistic mean fields @. A covariant treatment of 
nuclear matter provides a distinction between scalar and 
four- vector nucleon self energies, leading to a very natu- 
ral saturation mechanism. 

An example of a covariant EDF is known as the Rel- 
ativistic Mean Field (RMF) framework [3|. The corre- 
sponding effective Lagrangians provide a quantitative de- 
scription of a variety of ground-state data (masses, charge 
radii, ...). RMF, however, does not consider explicitly 
the Fock term. The exchange contributions are implicitly 
taken into account through the fit of model parameters to 
structure data. A more involved approach, the Relativis- 
tic Hartree-Fock (RHF) theory 0] , includes the exchange 
contributions explicitly. Early RHF models used to pre- 



dict nuclei that were considerably under-bound compared 
to experiment. The reason was the lack of a medium 
dependence in the corresponding effective nucleonic in- 
teraction [3]. The effect of the nuclear medium was first 
taken into account by adding self-interaction terms for 
the a meson field @. Although some improvement was 
obtained, the RHF results were still not on the level of 
RMF model predictions. An explicit nucleon-density de- 
pendence of the nucleon-meson couplings was included in 
Ref. [Ij]. This brought a significant improvement, so that 
current RHF models provide a quantitative description 
of nuclear properties with a similar accuracy as the stan- 
dard RMF approach JlQ|. In particular, recent studies 
by W.H. Long et. al. [sRll] and H. Liang et. al. [H OJ] 
have shown that, compared to the RMF approach, the 
explicit treatment of Fock terms can improve the descrip- 
tion of nuclear matter and finite nuclei. Moreover, it 
explicitly takes into account the tensor contributions to 
the inter-nucleon interaction generated by the exchange 
of the 7T and p mesons. These contributions have been 
found to play an important role in the description of the 
evolution of shell structures in the framework of the shell 
model jl 51 ] . The explicit treatment of exchange contribu- 
tions in covariant EDF models enables the inclusion of 
the pion field, which contributes only via its Fock term, 
and the tensor p-nucleon coupling, which contributes pre- 
dominantly in the exchange channel. The pion contribu- 
tion is expected to improve the predicted evolution of 
shell structure [ll[ , whereas the inclusion of the tensor p 
-nucleon coupling cures artificial shell gaps that arise in 
covariant EDF based model [l2| . 

Another benefit brought by the explicit treatment of 
the Fock term deals with the RPA description of collec- 
tive excitations. For example, the RHF+RPA approach 
provides a fully self-consistent description of charge- 
exchange excitations [lH, in contrast to the RMF+RPA 
model in which additional terms have to be intro- 
duced [3. 

When considering applications an important challenge 
for the framework of EDF is the systematic treatment 
of collective correlations related to restoration of broken 
symmetries and fluctuations in collective coordinates. A 
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static nuclear EDF is characterized by symmetry break- 
ing - translational, rotational, particle number, and can 
only provide an approximate description of bulk ground- 
state properties. To calculate excitation spectra and elec- 
tromagnetic transition rates in individual nuclei, it is nec- 
essary to extend the self-consistent mean-field scheme to 
include correlations that arise from symmetry restora- 
tion and fluctuations around the mean-field minimum. 
RMF-based models have recently been developed that 
include the explicit treatment of collective correlations, 
and employed in spectroscopic studies of a variety struc- 
ture phenomena related to shell evolution (l7j . 

The RHF framework has so far been limited to the 
description of spherical nuclei. In this work we con- 
sider an extension of this approach to deformed, axially- 
symmetric nuclei, and introduce the Relativistic Hartree- 
Fock-Bogoliubov model with density-dependent meson- 
nucleon couplings (RHFBz). In Sec. |TT]the general for- 
malism of the RHFBz model is presented. In Sec. IIIII we 
present and discuss the first applications of the RHFBz 
model to ground-state properties of carbon, neon and 
magnesium isotopes. To compare the results with those 
obtained with the standard RHB model Q, calculations 
are first performed without the inclusion of the pion field 
(PK02 parametrization [H,[l9j]). The effects induced by 
the pion field are analyzed in RHFBz calculations based 
on the PK03 parameter set (lallSJ- Finally, Sec. lIVI con- 
tains a short summary and discussion of possible future 
studies. 



Vectors in isospin space are denoted by arrows, and 
boldface symbols will indicate vectors in ordinary three- 
dimensional space. The Dirac spinor ip denotes the nu- 
cleon with mass M. m a , m u , m p , and m^ are the 
masses of the cr-meson, the w-meson, the p-meson and 
the 7r-meson, respectively. g a , g u , g p and f„ are the 
corresponding coupling constants for the mesons to the 
nucleon. e 2 /4ir = 1/137.036. The (density-dependent) 
coupling constants and meson masses are parameters, 
adjusted to reproduce nuclear matter properties and 
ground-state properties of finite nuclei. Vt^ v , , and 
J- 1 *" are the field tensors of the vector fields u, p, and 
of the photon Q . A nucleon- density dependence of the 
meson-nucleon couplings accounts for medium polarisa- 
tion and three-body correlations 0, HH • The effective 
Lagrangian is, therefore, characterized by 8 free param- 
eters: 

m a , g ,y{p S at), b ai d a 

9u(Psat)> b u> 
a p 

a* (2) 

that are adjusted in a fit to experimental masses of twelve 
spherical nuclei ( 16 0, 40 Ca, 48 Ca, 56 Ni, 68 Ni, 90 Zr, 116 Sn, 
132 Sn, 182 Pb, 194 Pb, 208 Pb, 214 Pb), as well as to nu- 
clear matter properties (saturation point, incompressibil- 
ity modulus ifooj and symmetry energy at saturation J). 

The single-nucleon Dirac equation is derived by varia- 
tion of the Lagrangian (JTJ with respect to ^ 



II. FORMALISM OF THE RHFBZ MODEL 



[£7% - M - E] i>{x) = , 



(3) 



A. Energy Density Functional 

1. Effective Lagrangian and equations of motion 

The RHFBz approach is based on a phenomenologi- 
cal Lagrangian density formulated in terms of relevant 
degrees of freedom for nuclear structure, namely nucle- 
ons and mesons. Nucleons are treated as point-like Dirac 
particles. The effective in-medium interaction between 
nucleons is described by meson exchange, whereas the 
Coulomb interaction between protons is taken into ac- 
count by the electromagnetic 4-potential . The La- 
grangian density reads : 



= |i7 M 9p - M - g<j{pv)v - guiPv)!^ 



where £ stands for the nucleon self-energy. When the 
variation is taken with respect to the boson fields, a set 
of inhomogeneous Klein-Gordon equations is obtained: 



(□ + ml) a = -g^ifj 
(□ + ml) = g u tfi> 
(□ + m 2 ) pf> = gptp-y^rip 

(□ + ml) 7r = ^d„ [^V^] 

771 tt 



(4) 
(5) 
(6) 

(7) 
(8) 



where the conservation of the baryonic current — 
■ip-f^if and the Coulomb gauge choice (d^A^ = 0) have 
been taken into account. The Hamiltonian of the model 
is derived from a Legendre transformation of the La- 
grangian: 

H = (fxy 7 ; [-i.V + M] il> 



-{d^d^-ml^)--{F, v Fn 



(1) 



+ - / d x [g a a + g^l^ + gpl^-r 

-\ — — 757.V7f.f-l- ejuA 11 — — - ]tp . 
m-jr 2 



(9) 
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2. Inclusion of the Fock term 

To explicitly include the exchange contributions, it is 
convenient to eliminate the mesonic degrees of freedom 
in Eq. ([§]) using the formal solution of the Klein- Gordon 
equations (Eq. Q): 

H = J d 3 xij H 7 .V + M] if; 

+ \ j d3xd4y $( x )$(y) T ™(x,y)D m (x,y) ^(s/)V>(z) 

(10) 

where a summation over the repeated index to — 
{er, ui, p 7 ir, A} is implied. D m {x, y) represents the propa- 
gator of the boson to, whereas Y m (x, y) corresponds to 2- 
body interaction matrices [7[. The nucleon field is quan- 
tized, and the no-sea approximation [22] is adopted for 
the nucleon states. The nucleon field operator -0 can 
be expanded on an auxiliary one-body operator basis 

i 

i 

The Hamiltonian (ITU1) consequently takes the form: 

H = T + Y J V m (13) 

m 

where 

T = J2 c i c i [ ^ HT V + M] fj 

V m = \Y1 C l C ] C k c l J d 3 Xid 3 X 2 ]i{xx)fj{x2) 

(14) 



i,j.k,l 

T m (l,2)D m (x 1 ,x 2 ) fk(x 2 )fi(xi) 



The effective inter-nucleon interaction contained in the 
Hamiltonian (| 13[) is designed to be used in the self- 
consistent mean-field approximation. The ground-state 
of the nuclear many-body system is, therefore, approxi- 
mated by a Slater determinant: 



*o>=n c !io) , 



(15) 



where | 0) represents the single-nucleon vacuum. The 
energy density functional is then obtained by taking the 
expectation value of the Hamiltonian (| 13[) in the ground- 
state Slater determinant (IT51): 



(16) 



£*** lp} = ($ \H[p] |$ 



In particular, the expectation value of the potential en- 
ergy operator generates a direct and an exchange term. 



This RHF functional can be written in terms of the one- 
body density operator represented by the matrix ele- 
ments Pij = ( $ | cj c * I $ o 



£ 



RHF 



[p, (j> m ) = Tr [(-17. V + M + T m <j> m ) p] 



+ \Tr [V [p] p^®p\±\J d 3 x [{d^ m f + ml] (17) 

Here r m represents the one-body vertex function Q • The 
trace operator involves a summation over space-time co- 
ordinates and Dirac indices. V is defined by its matrix 
elements Vijki- 

VijM = J d 3 xid 3 x 2 fi{xi)jj{x 2 )T m {l,2)D m {xi,x 2 ) 

fk(x2)fi(x!) (18) 
Finally, the tensor product corresponds to: 



(19) 



B. RHF equations for systems with axial symmetry 

The minimization of the energy functional (|17p with 
the constraint that the single-nucleon density matrix 
refers to a Slater determinant leads to the RHF equa- 
tions: 



[h[p},p]=0 

(- A +m 2 m ) <p m = ±Tr (T m p) , 



where 



h[p] 
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RHF 



5p 



(20) 



(21) 



In coordinate space the set of equations (1201 can be writ- 
ten as: 

{-ia.V + (3M*(r) + [V(r) + E R (r)] } Mr,q % ) + T t (r) 
= e i f i (r,q i ) (22) 

' (- A +ml) a(r) = -g a {p v ) p s (r) 
(- A+ml) u°(r) =g u (p v ) Pv(r) 
(-A+ml)p° 3 (r)=g p (p v ) p tv (r) ^ 
- A A°(r) = e p c (r) 

where q denotes isospin projection quantum number, and 
a = 7°7, ft = 7° are Dirac matrices. In the Dirac equa- 
tion (|2"2"|) : 

• M* = M + S(r) is the Dirac effective mass. 

• S(r) and V(r) are the Hartree terms, i.e. they rep- 
resent the direct contribution to the nucleon self- 
energy: 

S (r) = g a (p v ) a (r) (24) 
V (r) = gM cj° (r) + g p (p v ) p° 3 (r) r 3 + eA° (r) 

(25) 
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• S fi (r) denotes the rearrangement contribution. It 
can be divided into a direct term £jy(r) and an 
exchange term Sp(r). The direct contribution to 
the rearrangement term reads: 



(26) 



Taking the a meson as example, the exchange con- 
tribution to the rearrangement term reads: 

= ^6 qk , qi [^f k (qk)fl(qi)}{r) 
k,l Pv 



dr' { D a (r, r') (qi)fk (q k )} (/)}(27) 



• T i{r) denotes to the Fock terms, i.e. the exchange 
contribution to the nucleon self-energy. For in- 
stance, the Fock term associated to the a meson: 



[g.fMfiiqW^PbafA^qj) 



(28) 



The sources of the inhomogeneous Klein- Gordon equa- 
tions (1231) read : 



Ps(r) = ^2fi(r)fi(r) (29) 

i 

Pv(r) = J2fi(r)fi(r) (30) 

i 

Ptv(r) = Y,fHr)-nfi(r) = t?; oton (r) - p™ utron (r) 

i 

(31) 

Pc(r) = £ fKr) 1 —^!^) = P r ton (r) (32) 



(33) 



In the case of deformed nuclei characterized by axial sym- 
metry, the label i of the single-nucleon wave function 
fi{r) refers to the set of quantum numbers: 



i = (n,n, g ) 



(34) 



where f2 denotes the projection of the total angular mo- 
mentum on the symmetry axis, II is the parity, and q the 



isospin projection that distinguishes protons and neu- 
trons. In cylindrical coordinates (r±,4>,z) the nucleon 
wave function takes the form [231 : 



Mr) 



/2vr 



iG+(r ± ,z;q t )e^~^ 
\iGr(r ± ,z; qi )e^+^ J 



(35) 



The RHF equations are solved by expanding the nu- 
cleon spinors and meson fields in the basis of a deformed 
harmonic oscillator. The eigenfunctions $ a (r),a = 
{n z ,n r ,mi,m s } of the deformed harmonic oscillator po- 
tential 



z (r ± ,z) = hlojiri + ±Mu 2 z z 2 (36) 



are expressed in terms of Laguerre and Hermite poly- 
nomials. They form a basis on which the nucleon wave 
functions /;(r) are expanded [23| : 



Fi(r, qi ) 



1 



/27T 



F+(r ± ,z;q t )e 1 ^-^ 



= £if ) te)</>a(r±,z,^) 



(37) 



Gi(r,q. 



1 



/ j -JO. 

a 



'2tt 

9a'(qi)<t>&(r_L,z,(p) 



G+{r ± ,z; qi )e 1 ^-^ 
Gr(r^z- qi y^ + ^ 



(38) 



In a deformed harmonic oscillator basis, therefore, the 
solution of the Dirac equation (|22|) corresponds to a di- 
agonalization of the matrix: 



Ba. 



= e« 



f { :\q.) 

9%\qi) 



(39) 



Detailed expressions for the Fock contribution to the ma- 
trices A, B, and C are given in appendix YK\ 



C. Pairing correlations 

For a quantitative analysis of open-shell nuclei, both 
spherical and deformed, it is necessary to consider also 
pairing correlations. The nucleonic pairing is treated 
in the context of the Bogoliubov framework [24j. The 
resulting RHFB model provides a unified description of 
particle-hole (ph) and particle-particle (pp) correlations 
on a mean-field level by using two average potentials: 
the self-consistent mean field that encloses all the long 
range ph correlations, and a pairing field A which sums 
up the pp-correlations. Pairing correlations in nuclei are 
restricted to an energy window of a few MeV around the 



Fermi level, and their scale is well separated from the 
scale of binding energies, which are in the range of sev- 
eral hundred to thousand MeV. There is no empirical 
evidence for any relativistic effect in the nuclear pair- 
ing field A and, therefore, a hybrid RHFB model with 
a non-relativistic pairing interaction can be formulated. 
Similar to most applications of the RHB model the 
central part of the Gogny force [25| will be employed in 
the particle-particle (pp) channel: 



V pp (l,2) 



E 



-((i-i-r^Ax,) 2 



B,P a 



i=l,2 

-H^ - MiP°P T ) , 



(40) 



with the set DIS [26] for the parameters p,i, Wi, Bi, 
Hi, and Mi (i = 1,2). A basic advantage of the Gogny 
force is the finite range, which automatically guarantees 
a proper cut-off in momentum space. 



20 oscillator shells is used for the solution of the Klein- 
Gordon equations in the mesonic sector for both RHFB 
and RHB models. 
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III. RESULTS AND DISCUSSION 

The explicit treatment of exchange contributions re- 
quires the calculation of non-separable two-dimensional 
integrals I a pp' a ' in momentum space (cf. appendix IA1 . 
These integrals involve a boson propagator and the func- 
tions Q a ,p (Eq. (|A12jl ). The Q functions are non- 
separable two-dimensional spatial integrals. Their nu- 
merical evaluation imposes considerable constraints on 
the size of the basis of a deformed harmonic oscilla- 
tor. We have verified that a RHFBz calculation with 
6 fermionic shells yields reliable results for ground-state 
properties of nuclei up to Z = 30 (zinc isotopic chain). 
Various ground-state quantities (mass, axial deforma- 
tion parameter, charge radius, chemical potential, single- 
particle energies,...) obtained from an expansion of the 
nucleon wave functions in a basis of 6 oscillator shells, 
agree within 1% from those calculated with an expan- 
sion in 8 shells. Fig [1] displays the binding energy of 
the deformed 20 Ne nucleus with respect to the number of 
fermionic shells calculated with the PK02, DDME2 and 
Gogny DIS effective interactions. The N s h e a = 4 cal- 
culation is unphysical but emphasizes the similar quali- 
tative evolution of both PK02 and DDME2 binding en- 
ergies with the number of shell up to N s h e ii — 8. We 
estimate a 1% numerical error for the PK02 observables 
obtained from a 6 shells calculation, as seen on FigQ] We 
have also compared the results for 16 and 40 Ca calcu- 
lated with the expansion in 6 oscillator shells to those 
obtained by solving the spherical RHF equations in co- 
ordinate space, discretized on a mesh of R m ax — 20 fm 
and with a step a = 0.1 fm [12]. The corresponding 
ground-state quantities and single-particle energies dis- 
play relative variations of less than 1%, validating the 
choice of 6 fermionic shells for model calculation of light 
nuclei. The results shown hereafter correspond to calcu- 
lations with 6 fermionic shells for the RHFB model and 
12 fermionic shells for the RHB model. An expansion in 



FIG. 1: (Color online) Evolution of the calculated 
binding energy of the 20 Ne nucleus with the number of 
major shell. The theoretical values are obtained with 

the PK02 dHH, DDME2 i and Go g n y D1S 
effective interactions. 



27] 



A. Ground-state observables 

This section presents results of the first application 
of the RHFBz model in the calculation of ground- 
state properties of carbon, neon and magnesium nuclei. 
Masses, radii and shapes are fundamental characteristics 
of nuclei and their description presents a basic test of any 
model and effective force. Therefore, the RHFBz model 
with the PK02 and PK03 effective interactions pj, El] 
in the particle-hole channel, and the central part of the 
Gogny DIS force [26[ in the particle-particle channel, 
are used to calculate densities, masses, two-neutron drip- 
lines, deformations and charge radii of the Z = 6, 10, 12 
isotopic chains. The PK02 parametrization corresponds 
to a covariant EDF that does not include explicitly the 
pion field. Thus, the effects of the one-pion exchange 
is taken into account implicitly through the fit of the 
model parameters to data. PK03 parameterizes a co- 
variant EDF that explicitly includes the pion degree of 
freedom. Neither PK02 nor PK03 include the tensor p- 
nucleon coupling. The values of the PK02 and PK03 pa- 
rameters are listed in Table |H together with those of one 
of the most successful RMF functionals: DD-ME2 [Hj], 
that has extensively been used in applications of the RHB 
model. In this section, to make a first study of the influ- 
ence of the Fock term on a similar ground, we compare 
RHFBz results obtained without the inclusion of the pion 
field (PK02 effective interaction) with those of the axial 
RHB model (DD-ME2 effective interaction). 
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PK02 


PK03 


DD-ME2 


m CT (MeV) 


534.461792 


525.667664 


555.1238 


m u (MeV) 


783.000000 


783.000000 


783.0000 


m p (MeV) 


769.000000 


769.000000 


769.0000 


gcr(Psat) 


8.920597 


8.895635 


10.5396 




10.550553 


10.802690 


13.0189 


9p{Psat) 


2.163268 


2.030285 


3.6836 




0.000000 


1.000000 


0.0000 




1.375772 


1.244635 


1.3881 


b a 


2.064391 


1.566659 


1.0943 


Ccr 


3.052417 


2.074581 


1.7057 




0.330459 


0.400843 


0.4421 




1.451420 


1.245714 


1.3892 




3.574373 


1.645754 


0.9240 


Cu 


5.478373 


2.177077 


1.4620 


d w 


0.246668 


0.391293 


0.4775 


a p 


0.631605 


0.635336 


0.5647 




0.000000 


0.934122 


0.0000 



TABLE I: Parameters of the PK02, PK03 0,113 and 
DD-ME2 [H] effective interactions. 



Fig.[2]and Fig. [3] compare the proton and neutron den- 
sities of neon isotopes calculated with the PK02 and DD- 
ME2 effective interactions. It appears that both models 
predict rather similar shapes for nuclei with A < 26, 
whereas for heavier Ne isotopes larger deformations are 
calculated with the RHB model with DD-ME2, especially 
for proton densities. 

In Fig. [4] we display the absolute deviations of the cal- 
culated binding energies from the experimental values of 
neon isotopes for the two relativistic effective interactions 
PK02, DD-ME2, the Gogny force D1S, and the Skyrme 
interaction SLy4. Positive deviations correspond to un- 
derbound nuclei. One might notice that deformed RHFB 
calculations with PK02 predict binding energies with a 
level of agreement with data comparable to that of the 
Gogny D1S force and, for heavier isotopes, slightly better 
than RHB with DD-ME2. Much larger deviations from 
data are calculated with the Skyrme force SLy4. Similar 
results are found in the carbon and magnesium isotopic 
chains. This is quantified in Table ILTl where we compare 
the root-mean-square deviations of theoretical binding 
energies for the Z — 6, 10, 12 isotopic chains. PK02 pre- 
dictions are closest to the experimental values in the neon 
and magnesium isotopic chains, whereas DD-ME2 gives 
the smallest rms deviation for carbon nuclei. 

The two- neutron separation energy S^n = E to t(Z, N) — 
Ekot(Z,N-2) of Mg isotopes, calculated with PK02 and 
DD-ME2, are compared to data in Fig.[SJ In general, the 
RHFBz results obtained with the PK02 parameter set 
are closer to the experimental two-neutron separation en- 
ergies. The last two- neutron bound Mg nucleus is 38 Mg 
in the RHFBz calculation with PK02, whereas the RHB 
model with DD-ME2 predicts 40 Mg to be the last bound 





PK02 


DDME2 


Gogny D1S 


Skyrme SLy4 


o c (MeV) 


2.144 


1.443 


3.185 


2.874 


O Ne (MeV) 


2.263 


2.429 


2.750 


4.342 


a Mg (MeV) 


2.480 


2.582 


3.337 


3.269 



TABLE II: Root-mean-square deviations from 
experimental data of binding energies calculated with 
the PK02 Opl, DD-ME2 Q, Gogny D1S [H and 
Skyrme SLy4 [29[ effective interactions from for the 
carbon, neon and magnesium isotopic chains. 



isotope. Compared to DD-ME2, the PK02 two-neutron 
separation energies are also found closer to data in the 
carbon and neon isotopic chains. 

The evolution of the axial deformation parameter (3 
with mass number along the carbon isotopic chain is 
displayed in Fig. [5] No experimental results extracted 
from B(E2) measurements are presented insofar as it is 
not adequate to directly compare the static deformation 
parameter from the dynamical one extracted from the 
experiment, in such light even-even nuclei. PK02 and 
Skyrme SLy4 predict deformations that are systemati- 
cally smaller than those obtained with the DD-ME2 and 
Gogny D1S interactions, or with the Skyrme SGII effec- 
tive force. In particular, PK02 and SLy4 predict basi- 
cally spherical shapes between 10 C and 16 C whereas, ex- 
cept for 14 C, rather large ground-state deformations are 
calculated with the other three interactions. The case 
of 16 C is particularly interesting. Early experiments at 
Riken indicated an anomalously small B(E2) [32j, and a 
strong prolate deformation |33j. Therefore, the nucleus 
16 C was thought to be characterized by valence neutrons 
decoupled from a quasi-spherical core. Different models 
corroborated these results (see for instance [34[ )■ Recent 
measurements of B{E2\2\ — > + ) in 16 C gave a value 
that is more consistent with what is observed in nuclei 
with similar Nf Z [35|, |36| . The current result does not 
support the description of 16 C in terms of valence neu- 
trons decoupled from the spherical core. The RHFBz 
model calculation with the PK02 interaction, in partic- 
ular, predicts the neutron and proton axial deformations 
in 16 C: p n = 0.08 and (3 p = 0.06, respectively. 

The charge radii of neon isotopes, calculated with 
PK02 and DD-ME2, are shown in comparison with data 
in Fig. [7] In lighter Ne nuclei the theoretical values pre- 
dicted by the RHB model with DD-ME2 are in much bet- 
ter agreement with experiment, whereas for A > 26 both 
models yield similar charge radii. It should be noted that 
contrary to DD-ME2, PK02 does not include the charge 
radii in its fit to spherical nuclei. 

Structure models can also be compared by considering 
the corresponding single- nucleon spectra. In Ref. [13 the 
predictions of the spherical RHF model were tested in 
comparison to data for 16 and 40 Ca. In Fig. [8] we com- 
pare the proton Nilsson orbitals of 28 Mg, calculated with 
PK02 and DD-ME2. Although the ordering of Nilsson 
states is the same for both interactions, i.e. the defor- 
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FIG. 2: (Color online) Proton density in the Ne isotopic chain. The single-nucleon densities calculated with the PK02 (RHFBz) and DD-ME2 (RHB) 
effective interactions are plotted in the (Oxz) plan with x, z S [— 6/m, 6fm]. The color code denotes densities in the interval [0/m -3 , 0.09/to~ 3 ] 
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FIG. 3: (Color online) Same as Fig. [2] for neutron density in the interval [0/m 3 , 0.1/m 3 ] . 
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FIG. 4: (Color online) Absolute deviations of the 
calculated binding energies from the experimental values 
of the Ne isotopic chain. The theoretical values are 
obtained with the PK02 [IM1, DDME2 Q, Gogny 
D1S 27] and Skyrme SLy4 [29j effective interactions. 
The data are from Ref. [3(j]. For 32 Ne the calculated 
values are compared to the extrapolated binding energy. 
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FIG. 5: (Color online) Two-neutron separation energy 
in the magnesium isotopic chain. The relativistic 
mean-field results: RHFBz with PK02 [II El, and 
RHB with DD-ME2 [Hj], are compared to data 
(Audi-Wapstra [H). 
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FIG. 6: (Color online) Axial deformation parameter /3 
of C nuclei as function of the mass number. The 
calculated values correspond to the PK02 [3 El, 
DD-ME2 [H, Gogny D1S, Skyrme SLy4 and Skyrme 
SGII [3l| effective interactions. 
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FIG. 7: (Color online) Charge radii in the neon isotopic 
chain. The theoretical values calculated with 
PK02 [H El and DD-ME2 are compared to the 
ISOLDE experimental values [33]. The error bars do not 
include the atomic factor and therefore should be 
increased by 10%. 



B. The PK02 versus PK03 parametrization 



mation is rather small, in general the density of states 
around the Fermi level is larger when calculated with 
the RHFBz approach. This originates from the larger 
effective nucleon mass characterizing the PK02 interac- 
tion compared to the DDME2 one. Namely, the den- 
sity of states depends on the effective mass [38j . which 
is increased by the spatial non-locality of the mean-field 
potential (Fock terms). 



In the framework of shell-model calculations the ten- 
sor contribution, arising from pion exchange, has been 
found to play an important role in the description of the 
evolution of shell structures with proton/neutron num- 
ber [II. The effect of including the pion field in the 
RHBz model can be analyzed using the PK03 effective 
interaction [HI, El . In the relativistic mean- field frame- 
work the pseudo-vector coupling of the pion to the nu- 
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FIG. 8: (Color online) Comparison of the single-proton 
levels of 28 Mg, calculated with the PK02 pj, Qj| and 

DD-ME2 [28| effective interactions. The levels are 
labeled by Nilsson quantum numbers. The dashed-linc 
denotes the chemical potential. 




FIG. 9: (Color online) Same as in Fig. [4] but with the 
PK03 RHFBz calculation in addition. 



cleon generates part of the tensor contribution to the ef- 
fective inter-nucleon interaction, the remaining part be- 
ing induced by the tensor p-nucleon coupling. Here we 
consider the differences between the PK02 and PK03 
parametrizations on binding energies, ground-state axial 
deformation parameters, and charge radii. Contrary to 
the PK02 effective interaction, the PK03 effective in- 
teraction explicitly includes the pion contribution during 
the fit to data. The other parameters of the Lagrangian 
are also affected (Table HJ , meaning that the inclusion 
of the pseudo- vector n — N coupling alters how correla- 
tions beyond mean-field are implicitly taken into account 
through the effective meson couplings. 

Figure [S] displays the absolute deviations of the the- 
oretical binding energies from data for the sequence 
of neon isotopes. In addition to the results shown in 
Fig. 21 here we also include the deviations obtained in 
the RHFBz calculation with the PK03 interaction. The 
PK03 results for the 18 Ne, 26 Ne, 28 Ne and 30 Ne are on 
the same level of accuracy or better than those obtained 
with PK02, i.e. without the explicit inclusion of the pion 
field, whereas they show a less good agreement with data 
for the other Ne nuclei. 

The evolution of the axial deformation parameter /3 in 
the neon isotopic chain is illustrated in Fig. [10] In gen- 
eral, the deformation predicted by PK03 is larger than 
that calculated with PK02 and, therefore, closer to the 
results obtained with the DD-ME2 and Gogny D1S ef- 
fective interactions. PK03 predicts an oblate shape for 
24 Ne (quasi degenerate in energy with a prolate solution 
at j3 = 0.3) whereas a prolate ground-state shape for this 
nucleus is obtained with PK02, DDME2, Gogny D1S 
(quasi degenerate in energy with an oblate solution at 
j3 = —0.15) and Skyrme SLy4 interactions. Moreover, all 
these interactions, except PK03 that predicts a prolate 
ground-state, give no deformation for 26 Ne and 28 Ne. 




FIG. 10: (Color online) Same as in Figure O but for the 
neon isotopic chain, and with the PK03 RHFBz 
calculations and with the PK03 RHFBz results with 
fir(p) sets to in addition. 



Finally, in Fig.[TT]we illustrate the differences between 
the PK02 and PK03 parametrizations on the calcu- 
lated charge radii of neon isotopes. The results obtained 
with PK03 are shown in comparison with the ISOLDE 
data |37[, and with theoretical values predicted by the 
DD-ME2 (RHB) and PK02 (RHFBz) effective interac- 
tions (cf. also Fig. [7]) . One might notice that the explicit 
inclusion of the pion contribution leads to an enhance- 
ment of the calculated charge radii as compared with the 
values obtained with PK02, bringing them closer to the 
predictions of DD-ME2 and, for the lighter isotopes, in 
better agreement with data. 
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FIG. 11: (Color online) Same as in Fig. [71 but with the 
PK03 RHFBz calculation in addition. 



C. The effect of including the pion field 

In order to isolate the effect of the pseudo- vector ir — N 
coupling on the observables, we compare calculations 
based on the PK03 effective interaction where the pion 
contribution is switched on and off. Figure [12] displays 
the absolute deviations of the theoretical binding energies 
from data for the sequence of neon isotopes. In addition 
to the results shown in Fig. [9j here we also include the 
deviations obtained in the RHFBz calculation with the 
PK03 interaction where the pion contribution is switched 
off. Switching on the pion coupling constant brings rele- 
vant binding to the Neon isotopes. 




FIG. 12: (Color online) Same as in Fig. [9l but with the 
PK03 RHFBz calculation with f n (p) = in addition. 



term predict a spherical shape. The effect of the pion held 
on single-nucleon spectra is illustrated in Fig. [T3] and [TU 
where we display the proton and neutron single-particle 
levels in 26 Ne obtained in RHFBz calculations with the 
PK03 effective interaction where the pion coupling is 
switched on and off. PK03 with the pion coupling set 
to zero, yields a spherical ground-state shape and the 
Nilsson levels are degenerate, in contrast the degeneracy 
is lifted in the calculation performed with the complete 
PK03 interaction. Here one notices a clear signature of 
the effect of the pion on single-nucleon spectra and the 
corresponding evolution of shell structures. 
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FIG. 13: (Color online) Comparison of the single-proton 

levels of 26 Ne, calculated with the PK02 effective 
interaction (on the left), the PK03 effective interaction 
where the pion coupling is switched on (on the right) 
and off (in the middle). The dashed-line denotes the 
chemical potential. 
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FIG. 14: (Color online) Same as in Fig. [T3l but for the 
neutron single-particle levels. 
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In Fig. [10] the comparison between the PK03 curve 
and the PK03 one where the pion coupling is set to 
zero shows that the prolate shape of 26 ' 28 Ne is driven by 
the pion, whereas interactions without the explicit tensor 



In Fig. [TTJ we illustrate the effect of including the 
pion field on the calculated charge radii of neon isotopes. 
Switching on the pion coupling constant results in a dras- 
tic enhancement of the Neon isotopes charge radii. 
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IV. CONCLUSION 



butions read: 



We have developed the relativistic Hartree-Fock- 
Bogoliubov model for axially deformed nuclei (RHFBz). 
An effective Lagrangian with density-dependent meson- 
nucleon couplings is used in the particle-hole channel, and 
the central part of the Gogny force in the particle-particle 
channel. The RHFBz quasiparticle equations are solved 
by expansion in the basis of a deformed harmonic os- 
cillator potential. The numerical complexity brought by 
the explicit treatment of the Fock term within relativistic 
mean field theory limits, at present, the size of the oscilla- 
tor basis for the expansion of the nucleon wave functions. 
The current version of the model provides a reliable and 
numerically stable description of ground-state properties 
up to the zinc isotopic chain. Further numerical opti- 
mization is possible and work is in progress to extend the 
size of the deformed oscillator basis to 12 fermionic shells, 
allowing a description of medium-mass and heavy nuclei. 
In this work illustrative RHFBz calculations have been 
performed for carbon, neon and magnesium isotopes. Re- 
sults obtained with the RHF effective force PK02 have 
been compared to experimental masses and charge radii 
and, in addition, ground-state deformation and single- 
nucleon spectra have been shown in comparison with the 
predictions of one of the most successful RMF meson- 
exchange interactions: DD-ME2, as well as with the re- 
sults calculated with the non-relativistic Gogny D1S and 
Skyrme SLy4 interactions. The effect of explicitly in- 
cluding the pion field has been investigated for binding 
energies, deformation parameters, and charge radii. The 
addition of the tensor p-nucleon coupling will complete 
the model and thus enable studies of the role of tensor 
components of the effective inter-nucleon interaction in 
the evolution of shell structures in deformed nuclei. 
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LZ(V) = (2m z + n r - rj) L™» - 2 (n r + m,) 

(A3) 

Taking the a meson as example, the exchange contribu- 
tion to the A, B and C matrices is given by the following 
expressions: 
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Appendix A: Explicit expressions for the exchange 
contribution 



K, al (ft) = - E K .* fT fe) >< 

5" 5" F . . - 



Each matrix A, B, and C in Eq. contains the ki- S" ma {2.msp){2.ms^ l )I'[„s l - l 

netic, direct, exchange and rearrangement contributions. "f' a ' S/3 ° 

The explicit expressions for the kinetic and direct contri- (A5) 
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The matrices A, B and C are obtained by replacing the 
integrals IZuv'u' an ^ ^uvu'a 1 m tne expressions for the 
corresponding matrices A, B and C, with the integrals 
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(A10) 

The Fock terms are characterized by non-degenerate con- 
tributions to the positive and negative ilj blocks. The 
integrals can be written as: 



dr' D m (r, /) [g m 0y 07' ] (r'OAl 6) 
These integrals can be written in the form: 
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Taking again the case of the a meson field as exam- 
ple, the exchange contribution to the rearrangement term 
reads: 

fdg a 
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In the basis of a deformed oscillator, relation (|A13j) takes 
the form : 
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